Excitability in optical systems close to Z 2 -symmetry 
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Abstract 

We report theoretically and experimentally on excitability in semiconductor ring lasers in order to reveal a mechanism of ex- 
citability, general for systems close to Z2-symmetry. The global shapes of the invariant manifolds of a saddle in the vicinity of a 
homoclinic loop determine the origin of excitability and the features of the excitable pulses. We show how to experimentally make 

£j-^a semiconductor ring laser excitable by breaking the Z2-symmetry in a controlled way. The experiments confirm the theoretical 

^ predictions. 
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^ Today, excitability in systems outside equilibrium is a very 
O fertile and interdisciplinary research topic. Excitability has 
been observed in various fields including Chemistry, Physics, 
Biology or neural networks for computation tasks |[l][2l[3]S[5] 
• |6]0- In particular, excitability attracts a lot of attention in the 
^ field of optics |l9l|T0l[8l[ni[l3[l5l[ni[l2l due to its application 
•h asa way to generate well defined optical pulses and therefore 
starting a quest for optical excitable units. Lasers with saturable 
r*| absorber |9, 10 1, optically injected lasers fT2"l l81. lasers with op- 
Q_itical feedback fPJI or VCSELs with opto-electronic feedback 
lfl6l have all been proposed as optical excitable units, witness- 
__, ing both the high interest on optical excitability as well as the 

difficulty in achieving an ultimate design. 
OO The aim of this paper is two-fold: we introduce a simple, in- 
tegrable and scalable optical excitable unit based on a semicon- 
ductor ring laser (SRL) and we disclose in general the excitable 
properties of the wide class of nonlinear dynamical systems 
with weakly broken Z2-symmetry ifTTl [181 . SRLs are semi- 
conductor lasers whose cavity has a circular geometry that can 
sustain operation in two counterpropagating directions, namely 
clockwise (CW) and counter-clockwise (CCW) [19] [20). It 
was shown both theoretically GTIl and confirmed experimen- 
tally 11221 l23l that the dynamical regimes of a symmetric SRL 
^ are described by a two-dimensional asymptotic system of the 
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form: 



= fi (0,1ft) (ir = / 2 (0,#), 



(1) 



with the vector field /12 (9, if/) being invariant under the trans- 
formation 9 — > -9, 0- — » 27r — 0-, making SRLs an optical pro- 
totype of Z2-symmetric systems. 
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For the majority of optical systems, excitability takes place 
around a homoclinic bifurcation of a stable limit cycle lfT2l [8j 
[2 [TO] [TTJ. The presence of the stable limit cycle can easily 
be observed in these systems as it leads to Q-Switching oscil- 
lations. However, it is known that only unstable limit cycles 
undergo homoclinic bifurcation in SRLs [24], which are intrin- 
sically more challenging to address as they are not associated 
to observable dynamical regimes. 

In the first part of this paper, we present the theory for ex- 
citability in SRLs. Our approach is based on the investigation 
of the global topology of the invariant manifolds in systems 
with a slightly broken Z2-symmetry. In such a way we can pre- 
dict the onset of excitability near the homoclinic bifurcation of 
an unstable limit cycle as well as the properties of single and 
multiple excitable pulses. In the second part of this paper we 
show how to experimentally break the Z2-symmetry in a SRL 
in a controlled way. Excitability is revealed, whose properties 
confirm the topological predictions. 

Consider a SRL operating in single-transverse and single lon- 
gitudinal mode. Two directional modes, CW and CCW, can 
operate in the ring cavity with different intensities Pcw/ccw and 
phases (pcw/ccw- A linear coupling parameter with amplitude K 
and phase <p K is used to describe the transfer of power between 
CW and CCW. The two directional modes operate in antiphase, 
conserving the total power Pew + Pccw EOl ED [221 [23], such 
that an asymptotic two-dimensional set of equations can model 
the SRL operation: 

9 = 7sin#cos# + 2(1 - 5)cos ((f> k + uj) 

-(1 -sinf7)[(l -S) cos (cf> k +ifr) (2) 
+ (1 + S) cos (<f> k - if/)] 
cos 9\p = a J sin 9 cos 9 

-(1 +6) (1 -sin 9) sin (fa- ip) (3) 
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Figure 1: Phase space portraits of Eqs. |2}-j3| for different values of the bias 
current /. Black lines depict the stable manifold of 5 , while the red (grey) lines 
represent the unstable manifold of the S . In (a) the unstable limit cycle is shown 
in the red (grey) dashed line. The parameters are the following: 5 = 4.5%, 
$k = 1.5. (a) J = 0.659, (b) J = 0.691. 



+ (1 - 8) (I + sin 6) sin (<p k + if/) 



Here the angular variable 8 = 2arctan y/Pccw/Pcw — tt/2 
quantifies the power partitioning between the counterpropagat- 
ing modes. The phase difference tfr = 4>ccw - <Pcw is the second 
dynamical variable in the system. J is the rescaled bias current 
and a is the linewidth-enhancement factor. The Z2-symmetry 
of Eqs. is continuously broken by the introduction of 

an asymmetry AK in the mode-coupling. The dimensionless 
parameter 6 = AK/2K measures the relative magnitude of the 
symmetry breaking, in what follows we take 6 = 4.5% which 
corresponds to breaking the symmetry of the system in favor of 
the CW mode. In order to have a quantitative comparison with 
the experimental results shown below, we choose K = 0.1715 
ns~ l as the dimensionless time in Eqs. <[2]»-(|3j» is scaled with the 
physical backscattering strength K. We now discuss the phase 
space topology when varying the bias current J. 

The phase space of the SRL consists of two stable states CW 
and CCW whose basins of attraction are separated by the two 
branches of the stable manifold of a saddle S Il23l l24l . The 
main changes in the topology take place when the current J 
crosses a critical value 7/ !om which corresponds to a homoclinic 
bifurcation of an unstable limit cycle as shown in Fig.[Tfa)-(b). 
For / < Jhom [Fig.[TJa)], a small unstable limit cycle surrounds 
CCW and prevents the operation of the SRL in the CCW mode. 
Both branches of the unstable manifold of S connect S with 
CW, whereas one branch of the stable manifold of S connects 
S with the limit cycle and the other branch with an unstable 
node in (0,0) [not shown]. Excitability is possible in this sce- 
nario when the system -initially residing in CW- is driven by 
a fluctuation beyond the stable manifold of the saddle point 
S and subsequently relaxes back following the unstable man- 
ifold of S [15]. The radius of the unstable limit cycle rapidly 
grows with the current until it disappears in a homoclinic loop 
at J = Jhom- For J > Jhom, the basins of attraction of the CW 
and CCW modes are separated by the two branches of the sta- 
ble manifold of the saddle, a scenario that suggests bistability 
[Fig. [TJb)]. One would expect the system to leave the basin of 
attractor of the CCW state and relax to the opposite mode fol- 
lowing a branch of the unstable manifold of S . However, close 
to the homoclinic bifurcation, the global shape of the invariant 
manifolds is reminiscent of the homoclinic loop. As it is shown 
in Fig. [T|b), the unstable manifold of S then spirals very slowly 
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Figure 2: Numerical solutions of Eqs. j2j-j3) revealing excitable pulses in time 
domain (a,c,e) and their projection on the phase space (b,d,f). (a) Single pulse; 
n: pulse width, (c) Double pulse; T pp : intra-peak interval, (e) pulse with 
two maxima associated with the "metastability" of the CCW state; T2 : distance 
between the two maxima. The parameters are the same as in Fig.^b). 



to the CCW mode. The distance between consecutive laps of 
the unstable manifold decreases when J — > Jhom- Close to the 
homoclinic loop, the system relaxes very slowly to the CCW 
state, therefore making the CCW operation metastable: noise 
fluctuations will remove the system from the basin of attraction 
of CCW before CCW operation can be established. 

In the same way, close to the homoclinic bifurcation, the dis- 
tance between the branches of the stable manifold of the sad- 
dle becomes negligible when compared to the diffusion length- 
scale induced by the noise. The generation of an excitable pulse 
therefore corresponds to a noise activated crossing of two arbi- 
trarily close thresholds. An example of an excitable pulse in 
time domain is given in Fig. |2|a) and the corresponding phase 
space trajectory is shown in Fig. [2jb). Opposite to the scenario 
in Fig.[TJa), here the unstable manifold of S plays no role in the 
excitability. Excitability disappears for / too small or too large. 
When / is too large, the metastability of the CCW turns into 
full stability such that the system becomes bistable. For J too 
small, the system will exhibit alternating oscillations [24|. 

For both / < 7/ !0) „ and J > Jhom, we expect the intra-spike- 
interval (ISI) between consecutive excitable events to be dis- 
tributed in an exponential Arrhenius way as for a standard ac- 
tivation problem (25). The duration t\ of the pulse is instead 
a deterministic time scale corresponding to a revolution of the 
system in the phase space [see Fig. |2|b)]. The peak- width has 
been calculated from numerical integration of Eqs. <[2j-(|3j by 
averaging the Full-Width-Half-Maximum of a sample of 771 
peaks yielding t\ = 27.54ns. 

The shape and the features of the excitable pulses close to Z2- 
symmetry can also be explained by the topology of the invariant 




Figure 3: Schematic of the SRL and the experimental setup. The notation is as 
follows. CW/CCW: counterpropagating modes; I wg bias current on the waveg- 
uide; I r : bias current on the ring. 

manifolds. Consider the double pulses shown in Fig.[2jc): two 
well defined pulses are emitted by the SRL with an average 
peak-to-peak interval t pp = 49.22ns (averaged on a sample of 
286 double pulses). Double pulses can be explained as follows. 
The emission of the first pulse is due to a noise induced acti- 
vation; during the pulse, the system evolves deterministically 
towards the CW state. When the Z2-symmetry is weakly bro- 
ken and J is close to 7/, om , the deterministic evolution of the 
pulse brings the system in the vicinity of the stable manifold of 
S . A second excursion in the phase space, without residence in 
the CW mode is therefore possible as shown in Fig.|2|d). The 
peak-to-peak interval for such double pulses is determined by 
the duration of the deterministic rotation around the saddle in 
the phase space which is approximately twice the pulse-width. 
This topological insight is confirmed by the almost 2 : 1 rela- 
tion between t pp = 49.22ns and the pulse-width ti = 27.54ns. 
We remark that deterministic mechanisms for the generation of 
multiple pulses such as those in Ifl2l |8) require bifurcation sce- 
narios that cannot exist in two-dimensional systems and there- 
fore cannot be observed for Eqs. (|2j-([3j and in general for Z 2 - 
planar systems. 

Finally, excitable events occur that are characterized by two 
maxima separated by a minimum that does not reach the ground 
state such as the one shown in Fig.|2|e). The presence of such 
pulses is evidence of the metastability of the CCW-mode as 
discussed above: if the system enters the basin of attraction of 
CCW, it starts a rotation around CCW which results in the os- 
cillation of power. Due to the closeness to the homoclinic bifur- 
cation, during the rotation the system remains close to the stable 
manifold of S . The noise eventually drives the SRL back into 
the basin of attraction of the CW mode as shown in Fig. |2jf). 
The interval between the two peaks corresponds to a rotation 
around the CCW node and can be estimate as being 1 /2 of the 
period of a full rotation around the phase space. 

The experiments have been performed on an InP-based 
multiquantum-well SRL with a racetrack geometry and a free- 
spectral-range of 53.6 GHz. The device operates in a single- 
transverse, single-longitudinal mode at A — 1 .56fim. The power 
is extracted from the ring cavity by directly coupling the ring 
to a bus waveguide, and collected at the chip facets with a 
cleaved optical fibre as shown in Fig. [3] An electrical contact 
has been applied to the bus waveguide which can be indepen- 
dently pumped. The purpose of the contacting is twofold: it 
allows to counteract absorption and to amplify on-chip the sig- 
nal emitted by the ring; more interestingly the presence of a 
contact allow us to continuously break the symmetry of the de- 
vice in a controlled way by changing in an asymmetric way the 



I r [mA] 


ISI [pis] 


n [ns] 


t pp [ns] 


t 2 [ns] 


44.07 


23.3 


30 


53.0 


25.0 


44.62 


95.20 


26.7 


51.0 


26.6 


45.07 


110.7 


32.5 


49.5 


25.5 


45.38 


154.6 


33.9 


X 


28.0 



Table 1 : Measured peak properties for different values of the bias current /,. on 
the ring. The current on the waveguide for I p = 45.07mA is I wg = 14.45mA; 
for the other bias currents I wg = 14.0mA. The uncertainty is ±0.01mA for the 
bias current and ±4ns for the other quantities. For /,- = 45.38mA no double 
peaks were observed. 

relative strength 5 and phase A0t of the coupling between CW 
and CCW. Using the cleaved facet of the fibre as a mirror, we 
are able to reflect power from one mode (for instance CCW) 
back into the waveguide and finally to the counterpropagating 
mode in the ring. The amount of power that is coupled to the 
CW mode can then be controlled by tuning the current I wg on 
the waveguide, whereas its phase can be tuned by positioning 
the fibre facet with a piezoelectric controller. The chip con- 
taining the SRL is mounted on a copper mount and thermally 
controlled by a Peltier element which is stabilized at a temper- 
ature of 21.00°C with an accuracy of 0.01°C. We analyse the 
output power of the CCW mode using a fast photodiode con- 
nected to an oscilloscope with a sampling rate of 4.0ns. This 
sampling rate is slow enough for the purpose of recording long 
time series for statistical purposes, and it is fast enough to sam- 
ple accurately the individual excitable pulses. 

Due to the intrinsic symmetry of SRLs, the device is expected 
to operate with equal probability in either the CW or the CCW 
state. However, by breaking the symmetry using the piezo- 
controlled optical fibre and the bias-current on the bus waveg- 
uide, the operation in one of the two states (for instance CW) 
can be favored [see Fig. |3j. Time series of the power emitted 
in the CCW direction are shown in Fig.^a) for a bias-current 
on the ring of I r = 45.38mA and for a waveguide current of 
I wg = 14.0mA. It is clear from Fig. |4|a) that the ring operates 
most of the time in the CW mode. Short excitable pulses such 
as the one shown in Fig.Qb) are observed. The average width 
of these pulses is 33.9 ± 4ns, which agrees with the simulated 
value ri = 27.54ns. The ISI are distributed exponentially (see 
inset in Fig.|4jb)) with an average value of 154/is. We can there- 
fore deduce that the pulses have the characteristic signatures of 
excitability, being generated in noise-activation process across 
a threshold and having thereafter a deterministic evolution. We 
remark here that no residence in the CCW state was observed. 

The dependence of the IS I and pulse width T\ onI r is shown 
in Tab. [T| the IS I increases with I r , consistent with the picture 
of a noise activation process, whereas the pulse width remains 
almost constant. 

Among the pulses in the time series, double-pulse events are 
observed such as the one shown in Fig.[4|c). This kind of events 
consists of two well defined consecutive pulses separated by a 
drop of power to the ground state. The peak to peak interval 
T pp for different values of the bias current is shown in Tab. [T] 
Other pulse-shapes such as the one in Fig.Qd) have also been 
observed. They are characterized by two well-defined maxima 
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Figure 4: Experimental characterization of excitability in SRLs. (a) Time trace 
displaying excitability for /,- = 45.38mA (b) detail of (a) showing a single 
excitable pulse; (c) double pulse excitability for /,- = 44.62mA; (d) metastabil- 
ity for I r = 44.62mA. Inset of (b): distribution of the intra-spike intervals for 
I T = 45.38mA. 



separated by a shallow dip, and are a signature of the metasta- 
bility of the CCW mode as discussed above. The interval T2 
between the consecutive maxima is shown in Tab. Q] 
We observe from Tab.[T]that the ratio between the experimental 
values of t\ and r pp is in agreement with the 1 : 2 ratio pre- 
dicted by our topological arguments. In the same way, the ratio 
between T2 and r pp agrees well with the 1 : 2 ratio predicted by 
the theory. 

In conclusion, we have investigated excitability for generic 
planar systems close to Z2-symmetry and disclosed how the 
shape of the invariant manifolds lead to noise-activated pulses. 
An optical excitable unit based on this concept has been imple- 
mented using a multi quantum-well SRL with slightly asym- 
metric mode-coupling. Such unit can in principle be integrated 
on chip and does not require external optical injection or feed- 
back from an external cavity. We used a topological analy- 
sis to predict the features of three different types of excitable 
pulses, as well as quantitative relations between the relevant 
time-scales. The predictions of the theory have been confirmed 
by the experiments. We have numerically verified that the use 
of other non-Z2 perturbations leads to similar results. 
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